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Abstract We consider the Berkson model of logistic regression with Gaussian and homosce¬ 
dastic error in regressor. The measurement error variance can be either known or unknown. 
We deal with both functional and structural cases. Sufficient conditions for identifiability of 
regression coefficients are presented. 

Conditions for identifiability of the model are studied. In the case where the error variance 
is known, the regression parameters are identifiable if the distribution of the observed regres¬ 
sor is not concentrated at a single point. In the case where the error variance is not known, 
the regression parameters are identifiable if the distribution of the observed regressor is not 
concentrated at three (or less) points. 

The key analytic tools are relations between the smoothed logistic distribution function 
and its derivatives. 

Keywords Logistic regression, binary regression, errors in variables, Berkson model, 
regression calibration model 
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1 Introduction 

Statistical model. Consider logistic regression with Berkson-type error in the ex¬ 
planatory variable. One trial is distributed as follows. is the observed (or as¬ 
signed) suiTogate regressor. The true regressor is + [/„, where the error 

Un ^ NiO, T^) is independent of The response is a binary random variable 
and attains either 0 or 1 with 


p(y„=i I 


exp(/3o + PiXn) 

1 -f exp(/3o +' 
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We consider both functional model and structural model. In the functional one, 
are nonrandom variables, and in the structural one, are i.i.d., and therefore 
in the latter model, y„) are i.i.d. random triples. 

The couples F„), n = 1,..., iV, are observed. Vector 0 = (/3o, /3i)^ is a 

parameter of interest. 

The error variance can be either known or unknown, and we consider both 
cases. The conditions for identihability of the model (or of the parameter /?) are pre¬ 
sented. 

Overview. Berkson models of logistic regression and probit regression were set up 
in Burr [1]. For probit regression, it is shown that the introduction of Berkson-type 
error is equivalent to augmentation of regression parameters. As a consequence, the 
Berkson model of probit regression is identihable if is known and is not identihable 
if is not known. 

The identihability of the classical model was studied by Kiichenhoff [3]. He as¬ 
sumes that both the regressor and measurement error are normally distributed. Then 
univariate logistic regression is identihable (here can be unknown), and multiple 
logistic regression is not identihable. Our results can be proved similarly to [3] if 
we assume that the distribution of the surrogate regressor has an unbounded 
support. 

For classihcation of errors-in-variables regression models and various estimation 
methods, see the monograph by Carroll et al. [2]. 

Identihability of the statistical model can be used in the proof of consistency of the 
estimator. For known r^, the strong consistency of the maximum likelihood estimator 
is obtained by Shklyar [4]. But if is not known, the maximum likelihood estimator 
seems to be unstable (see discussion in [2] or [3]). 


2 Convolution of logistic function with normal density 


Consider the function 


Lo(a;,a2) = ^ ^ iV(0,a^), a: G R, > 0, 

that is, Lq{x, 0) = e“/(1 -f e“) and 

^ y_oo 1 + exp(a; - f) 


( 1 ) 


Denote the derivatives w.r.t. x 

Qk 

Lk{x,a'^) = -^Lo{x,a‘^). ( 2 ) 

Differentiation of L)^{x, a^) with respect to the second argument is described in 
Appendix A. 

The distribution of Yi given X”*”* is 

p[y, = 11 = E[p[r, = 11 \ 
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= E 


exp(/3o + PiXj) 

1 + exp(^o + PiXi) 



(3) 


since [/?o + PiX, \ ^ A^(/3o + PiX°'-p PIt^). 


3 Identifiability when is known 

Theorem 1. If in the functional model not all are equal, then the model is 
identifiable. 

Proof. Suppose that for two values of parameters P^^'^ = and = 

{pf, pf), P^^'^ f P^‘^\ the distributions of observations are equal. Then for all 

P;30)m = 1) = = 1), 

L^{pf + pfx°'°f [pffr^) = Lo{pf + pfxf'^f [pffr^)- 

However, by Lemma 4.1 from [4] the equation 

Lo{pf + Pfx, {Pffr'^) = LQ{pf + P^fx, {pff'T'^) 

has no more than one solution x. Hence, all Xare equal. □ 

By definition the degenerate distribution is the distribution concentrated at a sin¬ 
gle point. For the next theorem, see the proof of Theorem 5.1 in [4]. 

Theorem 2 ([4]). If in the structural model the distribution ofX°^^ is not degenerate, 
then the parameter P is identifiable. 

4 Identifiability when is unknown 

For fixed the function Lq{x, a^) is a bijection R —(0, 1). Hence, for fixed af 
and cr|, the relation 

-^0(2/, CTi) = Lo(a::, 0-2) ( 4 ) 

sets the bijection R. —>■ R; see Fig. 1. 

Lemma 3. For fixed erf > 0 and x| > 0, the sign of the second derivative of the 
implicit function (4) is 

= sign(x 2 - Xi) signjx). 

Proof. Differentiating (4), we get 

Li [y, CTi) dy = Li (x, X 2 ) dx; 

^ _ Li{x,al) 
dx Li(j/,x2)’ 
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Fig. 1. The plot to equation Lo{y, af) — Lo{x, ctI) for af < crl 


Then 


d^y _ L2{x,al)Li{y,al) - Li(a;, a|)£2(t/, cti)^ 

dx2 Li{y,alY 

_ L2{x,al)Li{y,alY - Li{x,a‘^fL2{y,al) 

^ / L2{x,al) _ L2{y,ai) \ ^ Li(x,a|)^ 
\Li{x,aiy Li(y,alY) Li{y,al) ' 


Thus, 




L2{x,al) 

Li(x,ct|)2 


L2{y,<Ji) \ 

Li{y,<xlY )' 


(5) 


Denote by y{z,a'^) the solution to the equation Lo{fj,,a'^) = z. Note that as 
L[j{x,aY is the cdf of a symmetric distribution, sign(Lo( 2 ;, cr^) — 0.5) = sign(a;). 
Therefore, sign(/i(z, ct^)) = sign(z — 0.5). Find the derivative 


/ L2{y{z,v),v) \ 
dv \Li{fj,{z,v),vY J 

for fixed 2 . By the implicit function theorem, 

dy{z,v) _ L 2 {fJ.{z,v),v) ^ 

dv 2 Li(/i(z, u), u) ’ 


d f L 2 {x,v) \ _ L 3 (x,v)Li{x,v) - 2 L 2 {x,vY 
dx\Li{x,vY) Li{x,vY ’ 

d / L 2 {x,v) \ _ L 4 {x,v)Li{x,v) — 2 L 2 {x,v)Ls{x,v) 
dv\Li{x,vY) 2Li{x,vY 

Then 

d / L 2 {y{z,v),v) \ _ L 2 L 3 L 1 — 2 L 2 L 4 L 1 — 2 L 2 L 3 

dv\Li{iJ,{z,v),vY) 2Li L\ 2L\ 
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2L\ 


where Lk are evaluated at the point (/i(z, u), v). By Lemma 10, 
d / L2{fiiz,v),v) 


Sign 


du \Li{fi{z,v),v)‘^ 


^ = sign(/r(z, u)) = sign( 2 ; — 0.5). 


The function v i—>■ ^fpiz’v) iJp monotone (it is increasing for z > 0.5 and 
decreasing for z < 0.5). For x and y satisfying (4), 

X = cr|) and y = y[z, aj) 

with z = Lo{y, af) = Lo{x, ct|); note that sign(z — 0.5) = sign(a:). Then 
L2{x,al) L2(y,al) 


sign 


^Li(a;,cr|)2 Li{y,aiy 
and with (5), we can obtain the desired equality 


= sign((T| - af) sign(a;), 



sign(^) = sign(cr| - al) sign(a;). 

□ 

Lemma 4. 

The equation 



Lo{Pf +pfx, crj) = Lo{pf +pfx, al) 

(6) 

has no more than three solutions, unless either 



and a\ = cr| 

(7) 

or 

^W=/3f)=0 and Lo{pf,al)=L^{pf,al). 

(8) 


In exceptional cases (7) and (8), equation (6) is an identity. 


Proof. The proof has the following idea: if a twice differentiable function y{x) satis¬ 
fies (4), then the plot of the function either is a straight line (if a\ = (t|) or intersects 
any straight line at no more than three points. 

Consider four cases. 

Case 1. a\ = (t|. Since the function Lq{z, ct^) is strictly increasing in z, Eq. (6) 
is equivalent to 

+ ffx = 13 f + ffx. 

Equation (6) has only one solution if jdf ^ pf', it is an identity if and 

it has no solutions if pf = pf but ^ pf. 

Case 2. pf = 0 and pf ^ 0. For any fixed cr^, the function z i—>■ Lq{z, ct^) is 
a bijection R —(0, 1). Denote the inverse function p.{Z, a^): Lo{z, a^) = Z if and 
only if z = y{Z,a^). Equation (6) has a unique solution 

y[Lo{pf - pf 

" ■ Pf 
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Case 3. = (3^'^ = 0. Neither side of (6) depends on X. Equation (6) becomes 

\ <Ji) = Equation (6) either holds for all x or does not hold for 

any x. 

Case 4. a\ y (j\ and y 0. Make a linear variable substitution: denote = 
/3 q^^ + Then Eq. (6) becomes 

^ • {Z 2 - = Lo(z 2 ,cr|). (9) 

V /?! / 

Dehne the function 21 (^ 2 ) from the equation 

-^^o(zi(-Z2)) ^i) = Lo(^Z2,a'l). 

The function zi (Z 2 ) : K —K. is implicitly dehned by Eq. (4): there the equality holds 
if and only if y = zi{x). Hence, the function zi(z 2 ) satisfies Lemma 3. Equation (9) 
is equivalent to 

0 ( 1 ) 

21 (^ 2 ) - Po^'’ -^ • {Z2 - pf'^) = 0. (10) 

Pi 

By Lemma 3, 



Then the derivative of the left-hand size of (10) 

^(zi(^2)-4'^-^- (^2-/3^)^ (11) 

is strictly monotone on both intervals (—cxj, 0] and [0, -l-c»), and hence (11) attains 
0 no more than at two points. Then the left-hand side of (10) has no more than three 
intervals of monotonicity, and Eq. (10) has no more than three solutions. Equation (6) 
has the same number of solutions. □ 

Theorem 5. If in the functional model there are four different then the param¬ 

eters P and Pit'^ are identifiable. 

Proof. Suppose that there are two sets of parameters {P^^\ and (P^^\ 

that for a given sample of the surrogate, the regressors {Xon, n = 1,..., iV} provide 
the same distribution of F„, n=l, ... ,N. Then for all n = 1,..., A^, 


— 1) — P^(2)^(T-(2))2(i^n — 1); 
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The equation 

has at least four solutions. Then by Lemma 4 either 


or 

and 

( 12 ) 

In the latter alternative, 

(M'Y (r('))' = (/3py = 0 and pf = pf 

since Lo{bo, 0) = is a strictly increasing function in &o- CH 

Theorem 6. If in the structural model the distribution of Xq is not concentrated at 
three {or less) points, then the parameters P and Pfr'^ are identifiable. 

Proof. Suppose that there are two sets of parameters {P^^\ and 

for which the same bivariate distribution of Yi) is obtained. The random vari¬ 
able P[yi = 1 I satisfies Eq. (3) almost surely for each set of parameters. 

Hence, the equality 

L^ipf+Pfxf^f {pf)\T^^^f)=Lo{pf+pfxt’^ {pff{T^^^f) 

holds almost surely. The rest of the proof is the same as in Theorem 5. □ 

A Differentiation ofI/fc(tc,cr^) 

Consider the sum of two independent random variables ^ = A -f where A has the 
logistic distribution 


P(A < x) 


exp(a;) 

1 -f exp(a:) ’ 


cc S K, 


and ^ ^ N{Q, a^). We allow cr^ = 0, and then ^ = 0 almost surely. 

The function Lo(a;, ct^) defined in (1) is the cdf of Q, and the function Li{x, cr^) 
defined in (2) is the pdf of C,. 

The partial derivatives of Lk{x, v) are 


d d 1 

—Lk{x,v) = Lk+i{x,v), —Lk{x,v) = -Lk+ 2 {x,v); 


see the proof in [4, Section 2]. The functions Lk{x,v) are infinitely differentiable and 
bounded on R x [0, -foo). 
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Since the distribution of C is symmetric, 

Lk{-x,a^) = fc > 1 , 

that is, Li{x, a^) and ^ 3 ( 0 ;, cr^) are even functions in x, and L2{x, ct^) and ^ 4 ( 2 ;, a^) 
are odd functions in x. 


B The key inequality 

The next lemma is similar to Lemma 2.1 in [4]. Hence, the proof is brief; see [4] for 
details. 

Lemma 7. Let ^ and 77 be two independent random variables, where ^ ^ N{0, 1). 
Denote C = ^ ^ Pc(^) pdfofC,. Then 

^(Inpf(z)) = I C= 2 :], 
where 7 x 3 [p | C=z] third conditional central moment, 

Pair] I C=2:] = E[(r7- E[77 I C=2])^ | C=z] ■ 


Proof. We have 


Pc(^) = - v) 



Ee 




Then 


P(iz) 


p'(;{z) _ E[(p-z)e-^(^-^)'] _ 

Pciz) ~ ~ 




E7y2g-Kz-’;)^ £Q-hi^-’iy - 

X {E[p^{p - z)e-5(^-'/)'] (Ee“5(^-'?)')2 
+ Ep^e-il^-'))' E[(r7 - 

- 2 E[p(p - Epe-5(^-'')" 

- 2 E[(p - 

+ 2(Epe-5(^-’')')^ E[(p - z)e-5(^-’')"]) 
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(13) 


If r] has a pdf, the conditional pdf of rj given Q=z is equal to 


Pr,\C=z{y) 


Pri{y)e 

£e-k{z-riy ’ 


Otherwise, we can use the conditional density of rj w.r.t. marginal density 

dcdf^|,^^^(y) _ 
dcdf^(i/) 

Anyway, the conditional moments of rj given <^=z are equal to 


E[p'= I C=z] 


£Q-i{z-vy 


(14) 


From (13) and (14) it follows that 

^(Inp^(z)) = E[p3 I C=2] - 3E[p2 I C=2] E[?7 | C=^] + 2(E[?7 | (= 2 ])^ 

= P3[n\C=z]. □ 


Corollary 8. Let ^ and rj be independent random variables such that ^ ^ <7^). 

Denote C = C + ^> denote the pdf of by p(^{z). Then 

1 

— (Inpc(z)) = I C=^]- 

Lemma 9. Assume that the distribution of a random variable X satisfies the follow¬ 
ing conditions: 

1) X has a continuously differentiable density px(x). 


2) X is unimodal in the following sense: there exists a mode M G R such that for 
all X S R, we have the equality sign(p^(x)) = sign(M — x). 

3) Whenever xi < M < X 2 andpx{xi) = px(x 2 ), then px{xi) > —px{x 2 )- 

4) E|X|3 < 00 . 


Then p 3 {X) ■= E{X -EX)^ > 0. 
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Fig. 2. To proof of Lemma 9, part 1). Sample px(a;) and definition of xi{z) and X 2 {z) 


Proof. 1 ) EX > M. Denote by xi{z) and X2{z) the solutions to the equation 
px{x) = z (see Fig. 2 ): 


xi{z) < M < X2{z) ifO < 0 < max(px); 
xi{z) = M = X2{z) if z = max(pjf); 
Px{xi{z)) = Px{x2{z)) = z if 0 < z < max(px)- 


Represent the expectation as a double integral and change the order of integration: 

/ OO 

(x — M)px{x) dx 

-OO 

= M+ jj (x —M)dxdz 


{{x,z) I 0<z<px(a:^)} 
^max(px) / rX 2 {z) 


= M + 


>0 


( pX2{z) \ 

/ (x — M) dx I dz 
Jx^{z) ) 


= M+ i^2{z) - Mf - (M - Xi(z))^ 


(15) 


For all X2 > M, by the implicit function theorem, 

because px(a;i(px(2:2))) = Px{x2) implies p'^{xi{px{x2))) > -p^(x2) > 0 . 

Note that Xi{px{M)) = M. By the Lagrange theorem, 


xi{pxix 2 )) = M + {x 2 - M) ■ ■^xi{pxix 3 )) 


X3—M-\-{x2 — M)9 


for some 6 S (0,1); 

Xi{px{x 2 )) > M - {X 2 - M) forx 2 >M; 

xi(z) > M — (x 2 (z) — M) for 0 < z < max{px); 
X 2 (z) — M > M — Xi(z) > 0; 
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Case (16) 


Case (17) 




0 


The sign of f(t) 

Fig. 3. To proof of Lemma 9, part 2) 


{X2{z) - MY ^ {M - Xi{z)Y 
2 ^ 2 ’ 

the last integrand in (15) is positive, and then (15) implies EX > M. 

2) Consider the function 

f{t) =Px(EX + f) -px(EX-f), 

which is odd and strictly decreasing on the interval [— (EX — M), EX — M] . There¬ 
fore, f{t) attains 0 only once on this interval, that is, at the point 0 (see Fig. 3). 

If f > EX — M (more generally, \t\ > EX — M) and f{t) = 0, then /'(f) = 
p'xi^X + t) + p'xi^X — t) > 0 by condition 3) of Lemma 9. Therefore, f(t) 
can attain 0 only once on (E X — M, -foo), and if it attains 0 (say, at a point ti > 
EX — M > 0), it is increasing in the neighborhood of ty. 

Hence, there may be two cases of sign changing of /(f) (Fig. 3). Either 

3fi>0Va:GK : sign(/(f)) = sign(f) sign(|f| — fi), (16) 

or 

VxGK : sign(/(f)) = — sign(f). 


3) We have 


/ OO 

{x — EX)px{x) dx 

-CXD 

/ OO 

fpx(EX -f f) df 

-OO 

nOO pOO 

= / fpx(EX-f f) df-I- / (-f)px(EX - f) df 

Jo Jo 


(17) 
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(18) 

where f{t) is defined in the second part of the proof. 

Note that the case (17) is impossible because otherwise the last integrand in (18) 
would be negative and thus the integral could not be equal to 0. 

4) Similarly to (18), 

E{X-EXf= / 

Jo 

Subtract times Eq. (18), where ti comes from (16): 

poo 

E{X-EXf= / t{t^-tl)f{t)dt. 

Jo 

The integrand is positive for f > 0,t y fi, and hence/r 3 [X] = E(7f—EX)^ >0. □ 

Lemma 10. For all x gM. and > 0, 



sign(L4(a;, cr^)Li(x, — SL^^x, a‘^)L2{x, a^)Li(^x, ct^) + 2L2{x, cr^)^) 
= sign(x). 


Lemma 11 is needed to prove Lemma 10. The notation F{y) and i/q is common 
for Lemmas 10 and 11. 

For fixed x > 0 and a^, consider the function 


F{y) = In 


ey \ 

{ey + 1 ) 2 / 


(y - xf 

2a2 


(19) 


Its derivative 


F'{y) = 1-2 


ey 

ey + 1 


y-x 


a 


2 


is strictly decreasing, and 

lim F\y) = +oo, lim F'{y) = —oo. 

y—¥ — oo y—y-\-oo 

Hence, F'{y) attains 0 at a unique point. Denote this point by yo, and then 

sig-n{F'{y)) = - sign(j/ - yo). (20) 

Lemma 11. For the function F{y) defined in (19), /or yo satisfying (20), and for 
yo and j /4 such that F'[yf) + F'{yf) = 0 and yo < 1 / 4 , we have the following 
inequalities: 

1) 2/3 < 2/0 < 2/4 and F'{yo) = -F'{yfj > 0. 

2) 2/3 + 2/4 > 0. 

3) F'fyo) < F'fyi) < 0. 

4) F{yo) > F{y4^). 
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Proof. 1) The inequality ys < yo < yi is a consequence of (20), and (20) implies 
F'iys) > 0. 

2) y 3 +y 4 > 0. For all y GR, 

0,7' 

F'{y) + F'{-y) = ->0. 


Since F'{y 3 ) + F'^-y^) > 0 and F'{y 3 ) + F'{yi) = 0, we have F'{-y 3 ) > F'{y 4 ), 
and then — 2/3 < 2/4 because the derivative F'{y) is decreasing. 

3) F''(jj 3 ) < F"{y 4 ) < 0. The second derivative 


F"{y) = 


-2ey 

{ev + 1)2 


1 


cr 


2 


is an even function strictly increasing on [0, +c») and attaining only negative values. 
The inequalities 2/3 < 2/4 and 2/3 + 2/4 > 0 can be rewritten as I 2 / 3 I < 2 / 4 . and then 


F"{y3) = F"{\y3\) < F"(2/4) < 0. 


4) F{x 3 ) > F{x 4 ). Consider the inverse function 


(F') \t), fsK. 


Its derivative is 


Then 


d 

dt 




F''{{F')-^it)) 


±{F{{F')-\t)) - F{{Fr\-t))) 


F'iiFTHt)) _ t 
F"{{F')-^t)) F"{{F')-^{t)y 

t ^ -t 

F"{{F')-yt))^ F"{{F')-y-t)y 


Apply already proven part 3) of Lemma 11 . If 2 > 0, then (F')“^(2) < {F')~^{—t) 
(because {F')~^{t) is adecreasingfunction)andF'((F"')“^(2))+F'((F"')“^(—2)) = 

2 — 2 = 0. Then by part 3) 

F"{{F')~\t)) < F"{{F')~\-t)) <0, 2 > 0. 

Hence, 

±{F{{F')-\t)) - f((F')-\-2))) >0, 2 > 0. 

Note that 

i^((F')''(0))-2^((F')-'(-0)) =0. 

By the Lagrange theorem, for 2 > 0, 


F{{Fr\t))-F{{Fr\-t)) = t-^fF{{F')-\h))-F{{Fr\-h))) > 0 , 

( 21 ) 

where the derivative is taken at some point 2 i S ( 0 , 2 ). 

Substituting 2 = ^'( 2 / 3 ) > 0 (then —2 = F'{y 4 )), we obtain ^"( 2 / 3 ) — F{y 4 ) > 0. 

□ 








144 


S. Shklyar 


Proof of Lemma 10. Case 1. x > 0 and > 0. Recall that for fixed a^, Li{x, a^) 
is the pdf of ?7 + ^, where rj and ^ are independent variables, P{rj < y) = and 
^ ^ 7V(0, tr^) (see Appendix A). By Corollary 8 , 

d^ 1 

^(InLi(a;,cr^)) = I vH=x], (22) 


but 




(lnLi(x,cr^)) 


L^L‘^ — 3L^L2Li 27^2 

if 


(23) 


where Lk are evaluated at the point (x, a^). Since Li(x, tr^) >0, we have to prove 
that y, 3 [ri \ ry+^=x] > 0. Therefore, we apply Lemma 9. 

The pdf of the conditional distribution of 77 given 77 + ^ = x is equal to 


1 (y-x)^ 

Pri\ri+i=x{y) = (y-x )2 • 2 -" 


The pdf Prj\rj+^=x{y) is continuously differentiable. The conditional distribution has 

_ (y-x )2 

a finite fcth moment because y e is bounded for any fc G N. Hence, conditions 

1) and 4) of Lemma 9 are satisfied. 

Evaluate 


l^Pv\v+i=^(y) = 111( (ey^ 1 ) 2 ) - ^ 2 ^ - ln(Ee = F^y) + C, 

where the function F{y) is defined in (19), and C = — ln(Eexp(—depends 
only on x and and does not depend on y. 

We check condition 2) of Lemma 9: 

Pr,|r,+J=x(y) = 

^p^\^+^^x{y) = F'{y)e^^y'>+^-, (24) 

sign(^p,,|^+5=x(2/)) = sign(F'(?/)) = - sign(y - yo), 

and condition 2) holds with M = ?/o, where yo is defined just above (20). 

Now check condition of 3) of Lemma 9. The proof is illustrated by Fig. 4. Assume 
that Pr,\r,+^^xiyi) = Pr,\rj+i=x{y 2 ) and T/i <yo < 2 / 2 -Then F(yi) = F{y 2 ). 

Denote 

yA=iF')-\-F'iy,)). 

Then F'{yi) + F'{yi) = F'{yi) - F'{yi) = 0, and by (20), as yi < yo, we have 
> 0, F'{yi) <0,yi> yo > yi. By Lemma 11, F{yi) > F{yi). 

Hence, ^( 7 / 2 ) = F{yi) > Fijji). Because the function F{y) is decreasing on 
( 77 o,+oo) (see (20)), we have 7/2 < 7 / 4 . Since the function F'iy) is decreasing. 
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Fig. 4. To proof of Lemma 10. Checking condition 3) of Lemma 9 


F'{y 2 ) > F'{y 4 ) = -F'{yi), which implies F'{yi) + F'(t/ 2 ) > 0. By (24) we 

havep;|^^^^^(yi) > 0. 

All the conditions of Lemma 9 are satisfied. By Lemma 9 , y^[r] \-q -\- ^ = x\ > 0, 
and by (22)-(23), 

Li{x,a"^)Li{x,a'^Y — ?)L^{x,a^)L2{x,a^)Li{x,a^)+2L2{x,a‘^) >0 (25) 
for all a; > 0 and cr^ > 0. 

Case 2. a; < 0 and > 0. The distribution of ry + ^ is symmetric. Hence, 
Li(a;, cr^) and L^{x,a^) are even functions in x, and L 2 {x,a^) and L^{x,a'^) are 
odd functions in x. Then 

Li{x,a^)Li{x,a'^Y ~ '^L^{x,a^)L 2 {x,a^)Li{x,a^) + 2^2(2;,cr^)^ 

is an odd function in x. It is equal to 0 for a; = 0, and it is negative for a; < 0 by 
Case 1; see (25). 

Case 3. cr^ = 0. The function Li{x, 0) is the pdf of the logistic distribution, and 
Lk+i{x, 0) is its fcth derivative: 

qX /1 — pXX 

= (TT^ijs' ^2(1,0)= ; 


Then 


ULj - 3L3L2T1 + 2Li = f!!il^(-2e“); 

sign(L4L^ - 3L3L2T1 + 2L2) = sign(a;), 

where Lk are evaluated at the point (a;,0). 

Lemma 10 is proven. 


□ 
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